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ABSTRACT 

The Transit Timing Variations (TTVs) technique provides a powerful tool to detect additional 
planets in transiting exoplanetary systems. In this paper we show how transiting planets with 
significant TTVs can be systematically missed, or cataloged as false positives, by current 
transit search algorithms, unless they are in multi-transit systems. If the period of the TTVs, 
Pttv, is longer than the time baseline of the observations and its amplitude, Attv, is larger 
than the timing precision limit of the data, transiting planet candidates are still detected, but 
with incorrect ephemerides. Therefore, they will be discarded during follow-up. When Pttv 
is shorter than the time baseline of the observations and Attv is sufficiently large, constant 
period search algorithms find an average period for the system, which results in altered tran- 
sit durations and depths in the folded light curves. Those candidates can get subsequently 
discarded as eclipsing binaries, grazing eclipses, or blends. Also, for large enough Attvs, 
the transits can get fully occulted by the photometric dispersion of the light curves. These 
detection biases could explain the observed statistical differences between the frequency of 
multiple systems among planets detected via other techniques and those detected via transits. 
We suggest that new transit search algorithms allowing for non-constant period planets should 
be implemented. 
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In a transiting exo-planetary system, the presence of other plan- 
ets can be inferred by the effect of their gravitational pull on the 
transiting planet. That effect can be detected as periodic transit tim- 
ing variations (TTVs), with amplitudes potentially as large as sev- 
eral tens of minutes (e.glMiralda-Escude, 2002; Agol et al. [2005 
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Holman & Murray| [2005f|Heyl & Gladman| [20071 |Veras & Ford 



These theoretical predictions have been confirmed by two re- 
cent results from the Kepler satellite mission ( [Boruckiet al.||2004| ). 
In the first of those results, the system Kepler- 9, composed of two 
transiting Saturn-size planets, presents TTVs with average ampli- 
tudes of 4 and 39 minutes per orbit ( |IIolman et al.||20i0l l. The other 
reported system, Kepler-11, composed of five 2.3 -13. SM^nnA tran- 
siting planets within 0.25AU from the star, and a sixth more mas- 
sive planet in an outer orbit, shows TTVs amplitudes of up to 30 
minutes jLissauer et al.|[20TT| ). |Stefren et al.| ( |2010| l have also pub- 
lished five Kepler systems with more than one transiting object, 
although those have not been confirmed as planets. 
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It is surprising that TTVs have been first detected in multi- 
ple transit planet systems, i.e. systems where several planets are 
co-planar within a few degrees, since those are geometrically less 
likely than non-coplanar multiple planet systems where only one 
planet transits. However, no definitive detection of single-transiting 
planets TTVs has been reported. Some authors have reported hints 
of possible TTVs, forWASP-3b l |Maciejewski et al.||2010^, WASP 



10b jMaciejewski et al.| [20TT]|, WASP-5b (|Fukui etal] [20TT]), 



and most recently HAT-P-13b i Nascimbeni et al. 201 1; Pal et 



||DTT1) 

poll 



), but none of those claims has yet been confirmed. Ford 



FaLl 
et aL] 



li also recently reported over 60 TTVs candidate systems in 
the Kepler mission data. 

The apparent scarcity of single-transit multi-planet systems 
was also noticed by Fabrycky ( 2009 ), who showed how one third 
of the non-transiting planets discovered via the radial velocity tech- 
nique, and with orbital periods similar to the transiting planets 
known at the time, had companions versus no detected compan- 
ions for any of the transiting planets. Two single-transit planetary 
systems with a second companion have been detected since then, 
CoR oT- 7b (|Leger et al.| [2509i [Queloz et a\] [2009l l and HAT-P- 
I3b ( [Bakos et aL 2009^ , but the statistics still disagree. 

In this work we explore potential observational and transit 
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search systematics that could introduce single-transit multi-planet 
systems detection biases. Section 2 summarizes how the transit 
search algorithms used by successful surveys work and their poten- 
tial detection biases. Section 3 describes our model simulation for 
a generic single-transiting system with TTVs and its main observa- 
tional signatures. Section 4 describes the model simulation results. 
We discuss our findings in section 5. 



2 TRANSIT SEARCH SURVEYS AND ALGORITHMS 

Since the discovery of the first transiting planet HD209458b a 

decade ago ( [Charbonneau et aL||2000| [Henry et^L||2000) , several 

surveys have successfully detected many other transiting systems. 

All those surveys combined, i.e. the ground- based projects OGLE 

(Udalski et al. l99T), HAT (Bakos et al. "^OOTl, TrES (Dunham] 

fet al, 2004, Alonso et al., 2004), XO (McCuUough et al. 20d5l(, 

I 1 1 — ^ , I — 1 1 1 

SuperWASP (PoUac co et al.[ 2006) , and the recent space missions 

CoRoT ( [Auvergne et al.||2009| l and Kepler l |Borucki et "ar|[2004| ), 

have already discovered over 100 transiting planets. In addition, 

the Kepler mission has recently announced 1235 transit candidates 

( [Borucki et al. 201 Ibj , of which 90% are expected to be true exo- 

planets I Morton & Johnson||2011| l. 

Transit planet candidates in each of those surveys are iden- 
tified via custom-made search algorithms adapted to the specific 
characteristics of each survey (i.e. noise systematics, sampling rate, 
etc ...). All the algorithms follow common approaches, which ei- 
ther look at the unfolded light curves in search for transit-like sig- 
nals, or fold the light curves over trial periods until a transit-like 
signal is found. Because of the typically sparse sampling of sin- 
gle orbit light curves, the second approach is the one ground-based 
surveys generally use. All the above ground-based surveys use the 
Box-fitting Least Square (BLS) algorithm ( |Kovacs et al.[|2002| or 
some modified version of it. The BLS algorithm consists on fold- 
ing the light curves over a series of periods, and then fitting the 
folded curves with narrow box-shaped functions of difi'erent depths 
and widths until finding the best least-square fit. OGLE, TrES and 
HAT use this approach. XO uses a modified form of the BLS algo- 
rithm which includes an adaptive-mesh search routine to increase 
the speed at which optimal periods are identified ( |McCullough| 
|etal.[|2005) . SuperWASP also uses a modified form of BLS which 
selects the strongest peaks of the BLS periodogram via Newton- 
Raphson refinement l |Protopapas et al.[|2005) . They also refine the 
box-shape function approach by 1) estimating the planetary radius 
using the physical properties of the star and analytical approxima- 
tions between the transit depth and the stellar and planetary radii, 
and 2) inspecting the folded light curves for photometric ellipsoidal 
variations characteristic of stellar binaries ( [Collier Cameron et al.| 
[2006| >. 

The CoRoT transit search algorithms have been developed 
by difi'erent participating teams and applied independently to the 
data. The algorithms, described in detail in'Moutou et al.'f2005^, 
[Regulo et aL[ ( [2007| l, [Carpano & FrTd lund (2008) and Ren ner et al. 
([2008), follow approaches such as correlations with sliding transit 



templates, box-shaped signal searches, wavelet transformations or 
BLS. Some of the transit search routines in these algorithms work 
with folded data, while others use the unfolded light curves. How- 
ever, all the algorithms, except for the ones by Renner et al.[p008) 
and |Regulo et al. 1 2007| ), limit their analysis to strictly periodic sig- 
nals. The [Renner et al.[p008| l algorithm includes some tolerance 
for period variations by allowing period error margins of 0. 1 days. 
The[Regulo et al.[([2007[) algorithm, based on wavelet transforma- 



tions, can detect period variations of about 2%, depending on the 
precision of the data (Regulo, priv. comm.). Both algorithms use 
unfolded light curves. 

Finally, Kepler uses a custom-made transiting planet search 
(TPS) pipeline, which starts by applying a wavelet-based, adap- 
tive matched filter to minimize the light curves noise, followed 
by a transit waveform (i.e. box-shape-like function) to search for 
individual transit events in the unfolded data. Those events are 
then folded over a wide range of periods until the period that best 
matches the full set of events is found (see [Jenkins[|2002] [Jenkins] 
[eTaTl [20T0a|b{ [Borucki eTaTI [20 lla[ ). TPS, which is optimized to 
detect Earth-size planets, folds the transit candidate events at strict 
periodicities, not allowing for very large TTVs. However, because 
of the continued monitoring of targets, the pipeline does flag sys- 
tems with significant TTVs (e.g. Kepler-9, Kepler-11). 



3 MODEL SIMULATIONS 

We have conducted model simulations to assess the effectiveness of 
current transit search algorithms for the case of transiting planets 
with significant TTVs. Our model system consists of a Jupiter-size 
planet orbiting a Sun-like star, with an orbital period Pp = 3.6235 
days and impact parameter b = 0. The model transit, generated us- 
ing jktebop ( [Southworth et al.[[2004a|b) , and a linear limb darken- 
ing of 0.61, has a depth of 1.3%, easily detectable by both ground 
and space surveys. Notice that we use as example case a model with 
parameters typical of known transiting hot Jupiters for the purpose 
of illustrating the potential TTV detection biases. The effect is still 
the same for planets with smaller radius or larger transit impact pa- 
rameter, except for the fact that transit detection limits will depend 
more strongly on the noise levels of the data (i.e. shallower transits 
get diluted into the noise at increasingly lower values of Ajtv, see 
§4). 

Using that model we simulate a series of transits adopting as 
variable parameters the length of the observing blocks, the observ- 
ing cadence, the photometric dispersion of the light curves, ai,i„ 
and the amplitude and period of the TTVs, Ajjv and Pttv ■ In our 
simulations Ajtv is defined as a fraction of the orbital period of the 
planet, or equivalently, a fraction of the orbital phase. For the TTVs 
we assume a simple periodic model of the form 



To + nPp + ArTvPp sin 



2;r 



-nP, 



(1) 



where Pttv is given in days, n is the number of elapsed orbits with 
respect to an arbitrary initial transit epoch, Tq, and r,„,y is the mid- 
transit time of each transit. 

For the simulations we adopt values of Attv between 0.0025 
and 0.02, and values of Pttv between 25 and 360 days. In addition, 
to resemble as close to real observing conditions as possible, we 
perform the simulations with parameters typical of both ground- 
based and space-borne surveys. In the case of ground-based sur- 
veys we simulate two scenarios: the first one adopts an observing 
block length of about 3 months (90 days), a photometric disper- 
sion a pi, = 8 X 10"'', and an observing cadence of 8 minutes, which 
approximate the average values of those parameters for the ground- 
based surveys listed in section 2, over one observing season (year). 
In the second scenario we generate three observing blocks of 90 
days each over a three-year period to simulate long-term ground- 
based observations spanning over more than one observing season. 
In the simulations we also include the effect of day-night cycles, 
and a 70% probability of clear nights to estimate the total number 
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Figure 1. Effect of TTVs for 90-day observing blocks, when Pttv = 25, 100 and 360 days (columns left to right) and Ajtv = 0.0025, 0.01 and 0.02 (rows 
top to bottom). In each panel, the top curves are model data for cr,,;, = 10"'' (Kepler case), and the bottom curves are for o",,;, = 8 X 10"^ (ground-based case). 
All curves are folded over the best-fit periods. Solid circles show ground-based model data binned over 0.002 phase intervals. 



of observable transits. To simulate the conditions of a space-borne 
survey, we adopt parameters similar to the Kepler mission, i.e. a 
photometric dispersion cr,,;, = 10"*, a 30 minute cadence, and two 
different observing block lengths: one of 90 days, to compare with 
similar ground-based data, and another one with continuous obser- 
vations over 3 years, which is the nominal lifetime of the mission 
( |Borucki et ari |2004li. 

Each synthetic light curve is then folded over a wide range of 
trial periods, in a manner similar to the BLS algorithms, until the 
best period is found. That best period is designated as the one that 
minimizes the transit duration in the folded light curves, or equiv- 
alently, period that minimizes the transit's photometric dispersion. 
To interpret our results we use the raw space-borne model data, and 
the ground-based model data averaged over 0.002 phase bins. 



4 SIMULATION RESULTS 

The outcome of the simulations is illustrated in Figures 1 and 2. 
Figure 1 shows the results for the 90-day observing blocks, both 
for the ground-based and Kepler survey cases. Figure 2 shows the 
results for the 3-year observing blocks (3 seasons for ground-based 
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Figure 2. Same as in figure 1, but for the 3-year observing blocks. All plots 
correspond to Pttv = 100 days. The amplitudes of the TTVs represented 
are Attv = 0.0025, 0.01, 0.015 and 0.02. 



surveys and continuous observations for Kepler). In all cases, a 
BLS-type algorithm should detect transit-like dips for Attv < 0.02. 
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However, the characteristics of the curves change with the values 
of Attv and Pob^ i.e. the length of the observing blocks. 

Figure 1 reveals two distinct cases. When Pttv < Fobs (left 
column in figure), the full amplitude of the TTV is observed and 
any best-fit algorithm should fold the data about the correct orbital 
period of the planet. Pp. However, the transits widen and get shal- 
lower as A^jy increases. While the effect of these TTVs will be dis- 
cernible in Kepler-like data, those transiting planet candidates will 
be most likely labeled as false positives given the lower light curve 
data quality of a typical ground-based survey (signals that wide and 
shallow will be mistaken for grazing binary eclipses). Further tests 
show that for Attv s> 0.02 the transits get completely diluted in 
typical ground-based data. When Pjtv > Pohs (right column in fig- 
ure), the algorithms find a perfect transit shape, but with the data 
folded over an erroneous period 



P„\l+ InAr 



(2) 



which results from eq. 1 when Iny-^ ^ 1- This period offset also 
affects the Tq zero epoch value, and therefore the resulting transit 
ephemerides equation will be incorrect. Any follow-up observation 
will therefore miss the transit event and the candidate will be dis- 
carded as a spurious detection. 

For intermediate cases, where Pttv ~ Pots (middle column in 
figure), a best-fit folded light curve simultaneously yields a wrong 
ephemeris and shape-deformed minima, which may again elude 
correct detection. 

Figure 2 illustrates a situation similar to the left-side column 
of figure 1 (i.e. Pttv < Pobs), but for the 3-year observing blocks. 
This figure shows more clearly how the depth and width of the 
transits get altered as a function of Attv- In the case of Attv = 
0.01, for example, the folded ground-based light curve resembles a 
"V" shape, typical of binary star false positives. 



5 DISCUSSION 



[Terquem & Papaloizou] \2QQ1) and |Lo Curto et al.| J20T0| have 
suggested that the absence of multiple planet systems (and there- 
fore TTVs) among transiting planets discovered by surveys (in 
their majority hot Jupiters discovered from the ground until Ke- 
pler was launched) could be due to migration mechanisms. Mas- 
sive short period planets would have moved inwards disturbing the 
orbits of all other planets in their way. However, the recent detec- 
tion with Kepler of over 130 candidate systems showing some ev- 
idence of TTVs challenges those hypotheses ( [Ford et aL| |2011| ). 
Figure 3 represents the frequency of TTV candidates reported by 
|Ford et aLl ( |2011) in two separate histograms; one for planets with 
Rp < 6REanii, i-e. Neptune-size and smaller planets, and the other 
for planets with dREanh < Rp < iSREani,, i.e. gas giants. Given 
photometric precision and observing cycle limitations, gas giants 
with orbital periods of six days or less are the planets most likely 
to be detected from the ground. The bottom diagram in the fig- 
ure shows five TTV candidates that match those criteria. Assuming 
all of the [Ford et al.|J2011^ candidates truly are planets with TTV 
variations, and comparing those numbers with the total number of 
hot Jupiters candidates with P < 6 days in the current Kepler mis- 
sion's candidate list (~ 65 planets), 1 out of each 15 transiting hot 
Jupiters should show TTVs. Based on this estimate and given the 
number of hot Jupiter transiting planets currently confirmed, about 
7 systems with significant TTVs should have been already discov- 
ered from the ground. Doing a similar estimation with the num- 
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Figure 3. Histograms depicting the frequency of planets with TTVs in the 
jFordet al.H2011) sample as a function of orbital period. The top panel con- 
tains all planets smaller than about Neptune size (Rp < dREmih)', the bottom 
panel contains all Saturn-Jupiter size planets (6R£„,,;, <Rp < 15R£„,,;,.The 
vertical dashed line around P = 6 days indicates the, by definition, period 
limit of hot Jupiters. Most of the TTV detections for close-in planets cor- 
respond to the smaller planets, but there are also several massive close-in 
planets in the sample, in contrast with theoretical expectations. 



ber of multiple planet systems found to date via the radial velocity 
technique (source http://exoplanet.eu I, we find that 14 have planets 
with orbital periods P < 6 days and Msini > 0.5Mj. Of those 2 
are Jupiter-mass planets (assuming Msini > O.SMj), which gives 
a fraction of close-in massive planets in multiple systems of about 
1/7, although not all will present significant TTVs . Assuming also 
that a fraction of the |Ford et al.H20lT} TTV candidates will be false 
detections, we conclude that the true number of hot Jupiters show- 
ing TTVs falls somewhere below 1/15, and those systems are most 
likely being missed by current planet search algorithms. 

Our simulations reveal two effects that may explain why tran- 
siting planets with significant TTVs elude detection. One is finding 
the wrong transit ephemeris when Pttv > Pots- This effect does not 
preclude those systems from being labeled as transit candidates, but 
if a dubious candidate around, for instance, a relatively dim star is 
detected by a small telescope with a pi, > 0.01, and cannot be con- 
firmed by observing it at a higher photometric precision because 
the ephemerides are incorrect, the candidate will be discarded be- 
fore spectroscopic confirmation follow-up. 

The other effect, i.e. wider, shallower and "V"-shaped transits, 
may explain why planets with Pttv < Pots and significantly large 
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Figure 4. Representation of the transit deptli versus transit length for hypo- 
thetical planets with radius between 1 and 2 S^,,,, (solid line). The dashed 
line shows the same relation derived from the result of our simulations 
for Attv = 0.0025, 0.005, 0.01, 0.015 and 0.02, a 3-year observing block 
length, and Pjtv 100 days. 



Attv are lost. Figure 4 shows how, for our simulated \Rj„p planet, 
the depth and length of the transit changes in the folded light curve 
when Attv varies from 0.0025 to 0.02. The almost vertical line 
in the figure indicates, for comparison, how the transit length and 
depth of a real planet would vary when its radius increased from 1 
to 2 Rjup- Clearly, a modified BLS algorithm , or visual inspection, 
will find the depth and length of the candidate to be incompatible 
with the expected parameters of a real planet, and the system will 
be rejected as a grazing eclipsing binary. In the case of very low 
light curve scatter and continuous space-borne observations, like 
Kepler's, TTVs will be more difficult to miss, as show Figures 1 
and 2. However, faint targets will encounter the same problems as 
grotmd-based observations. Algorithms like the ones developed by 
[Regulo et aLl l |2007^ and |Renner et al.| ( |2008^ should have been able 
to detect TTVs in the CoRot data, but the limited monitoring per 
field (~ 150 days, |Auvergne et al.[ |2009l, and strong data system- 
atics may have prevented such detections. 

We conclude that, if our hypotheses are correct, the imple- 
mentation of new transit search algorithms optimized for the detec- 
tion of TTVs, and their application to existing ground-based transit- 
search databases, should lead to the detection of new single-transit 
systems with significant TTVs. 



ACKNOWLEDGMENTS 

We thank I. Ribas for useful discussions during the development of 
this work. We also thank J. Jenkins, D. Latham, C. Regulo and S. 
Renner for helpful clarifications about transit search algorithms of 
different surveys. 



References 

Agol E., Stefi'en J., Sari R., Clarkson W., 2005, MNRAS, 359, 567 
Alonso R., et al. 2004, ApJL, 613, L153 
Auvergne M., et al. 2009, A&A, 506, 411 



Bakos G., et al. 2009, in lAU Symp. Vol. 253 of lAU Symposium, 

HAT-South: A Global Network of Southern Hemisphere Auto- 
mated Telescopes, pp 354-357 
Bakos G. A., Lazar J., Papp I., Sari R, Green E. M., 2002, R\SR 

1 14, 974 
Borucki W., et al. 2004, in F. Favata, S. Aigrain, & A. Wilson ed.. 

Stellar Structure and Habitable Planet Finding Vol. 538 of ESA 

Special Publication, The Kepler mission: a technical overview. 

pp 177-182 
Borucki W. J., et al. 201 la, ApJ, 728, 1 17 
Borucki W. J., et al. 201 lb, ArXiv e-prints 
Carpano S., Fridlund M., 2008, A&A, 485, 607 
Charbonneau D., Brown T. M., Latham D. W., Mayor M., 2000, 

ApJL, 529, L45 
Collier Cameron A., et al. 2006, MNRAS, 373, 799 
Dunham E. W., Mandushev G. I., Taylor B. W., Oetiker B., 2004, 

PASP 116, 1072 
Fabrycky D. C, 2009, in lAU Symp. Vol. 253 of lAU Symp., 

What to Expect from Transiting Multiplanet Systems, pp 173- 

179 
Ford E. B., et al. 201 1, ArXiv e-prints 
Fukui A., et al. 201 1, PASJ, 63, 287 
Henry G. W., Marcy G W., Butler R. P, Vogt S. S., 2000, ApJL, 

529, L41 
Heyl J. S., Gladman B. J., 2007, MNRAS, 377, 1511 
Holman M. J., Murray N. W., 2005, Science, 307, 1288 
Holman M. J., others 2010, Science, 330, 51 
Jenkins J. M., 2002, ApJ, 575, 493 
Jenkins J. M., et al. 2010a, ApJL, 713, L87 
Jenkins J. M., et al. 2010b, in Society of Photo-Optical Instrumen- 
tation Engineers (SPIE) Conference Series Vol. 7740, Transiting 

planet search in the Kepler pipeline 
Kovacs G., Zucker S., Mazeh T, 2002, A&A, 391, 369 
Leger A., et al. 2009, A&A, 506, 287 
Lissauer J. J., et al. 201 1, Nature, 470, 53 
Lo Curto G, et al. 2010, A&A, 512, A48+ 
Maciejewski G, et al. 2010, MNRAS, 407, 2625 
Maciejewski G, et al. 201 1, MNRAS, 41 1, 1204 
McCullough P R., Stys J. E., Valenti J. A., Fleming S. W., Janes 

K. A., Heasley J. N., 2005, PASP 117, 783 
Miralda-Escude J., 2002, ApJ, 564, 1019 
Morton T. D., Johnson J. A., 201 1, ArXiv e-prints 
Moutou C, et al. 2005, A&A, 437, 355 
Nascimbeni V, Piotto G, Bedin L. R., Damasso M., Malavolta 

L., Borsato L., 201 1, ArXiv e-prints 
Pal A., et al. 201 1, MNRAS, 413, L43 
PoUacco D. L., et al. 2006, PASP 118, 1407 
Protopapas P, Jimenez R., Alcock C, 2005, MNRAS, 362, 460 
Queloz D., et al. 2009, A&A, 506, 303 
Regulo C, Almenara J. M., Alonso R., Deeg H., Roca Cortes T, 

2007, A&A, 467, 1345 
Renner S., Rauer H., Erikson A., Hedelt P, Kabath P, Titz R., 

Voss H., 2008, A&A, 492, 617 
Southworth J., Maxted P F L., Smalley B., 2004a, MNRAS, 351, 

1277 
Southworth J., Zucker S., Maxted P F L., Smalley B., 2004b, 

MNRAS, 355, 986 
Stefi^en J. H., et al. 2010, ApJ, 725, 1226 
Terquem C, Papaloizou J. C. B., 2007, ApJ, 654, 1110 
Udalski A., Szymanski M., Kaluzny J., Kubiak M., Mateo M., 

1992, AcA, 42, 253 
Veras D., Ford E. B., 2010, ApJ, 715, 803 



